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ABSTRACT 

We estimate the evaporation timescale for spherical H I clouds consisting of the cold neutral medium sur- 
rounded by the warm neutral medium. We focus on clouds smaller than Ipc, which corresponds to tiny H I 
clouds recently discovered by Braun & Kanekar and Stanimirovic & Heiles. By performing one-dimensional 
spherically symmetric numerical simulations of the two-phase interstellar medium (ISM), we derive the 
timescales as a function of the cloud size and of pressure of the ambient warm medium. We find that the 
evaporation timescale of the clouds of 0.01 pc is about IMyr with standard ISM pressure, p/ks ^ lO"^'^ K 
cm^-^, and for clouds larger than about 0. 1 pc it depends strongly on the pressure. In high pressure cases, there 
exists a critical radius for clouds growing as a function of pressure, but the minimum critical size is ^ 0.03 
pc for a standard environment. If tiny H I clouds exist ubiquitously, our analysis suggests two implications: 
tiny H I clouds are formed continuously with the timescale of IMyr, or the ambient pressure around the clouds 
is much higher than the standard ISM pressure. We also find that the results agree well with those obtained 
by assuming quasi-steady state evolution. The cloud-size dependence of the timescale is well explained by 
an analytic approximate formula derived by Nagashima, Koyama & Inutsuka. We also compare it with the 
evaporation rate given by McKee & Cowie. 

Subject headings: hydrodynamics - ISM: clouds - ISM: kinematics and dynamics 



1. INTRODUCTION 

Recently discovered "tiny H I clouds" have very small 
sizes, ~ 10^^ pc, and small H I column density, ~ 
lO^cm"^ jBraun & Kanekar 2005; Stanimirovic & Heiles 
This new population of clouds is providing a 
challenge to our understanding of the evolutionary cy- 
cle of the interstellar medium (ISM). The origin of the 
clouds is still under debate, but it is presumably the 
thermal instabihty (Field 1965; Balbus 1 986 )^n turbu- 
lent gas (e.g. VazQuez-Semadeni, Passot & Pououetl 11995 
|&)vama & Inutsuka 2002, 2004; Ki-itsuk & Norman 2002aT 
l^nnebelle & Audit 2005) and/or the fragmentation of cold 
clouds crushed by interstellar shocks (Nakamura et al. 2005|. 
Apart from the formation process, we can extract many in- 
formation about the surrounding ISM via investigation of the 
evolution and statistics of the clouds. In this Letter, we thus 
study the evaporation rate and timescale of the tiny H I clouds, 
which are key quantities characterizing the fate of the clouds. 

One of the simplest model is the clouds that con- 
sist of the cold neutral medium (CNM) with tempera- 
ture T ~ 10^~^ K, surrounded by the warm neutral 
medium (WNM) with T ^ lO'* K, under pressure equi- 
librium. The two phases are thermally stable balancing 
the heating rate due to, e.g., photo-electric heating of dust 
grains, with the cooling rate due to, e.g., line cooling 
by fine-structure transition of C II atoms. This is based 
on a widely accepted two-phase description of the ISM 
(field. Goldsmith & Habing 1969; Wolfire et al. 2003). The 
evolution of clouds is, therefore, described as the motion 
of the interface, or front, between the CNM and WNM, 
driven by the thermal conduction . Based on this p icture, 
IZeFdovich & Pikel'ned (fT969l) and lPenston & Brownl (fT97nh 
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considered isobaric, steady motion of fronts in plane-parallel 
geometry. They clarified that the motion is determined by 
pressure and that there is a saturation pressure for which a 
static front can exist. Their work has been substantially ex- 
tended by many authors (e.g. Elchick. Reaev & Sniegel 19911 
Ebhick. Reeev & Shaviv 1992£' (Ferrara & Shchekinov«ll993[ 
.Hennebelle & Perauh 1999) . 

Compared to the analysis of the frontal motion in plane- 
parallel geometry, the evolution of spherical clouds has been 
much less analyzed. Graham & Langer ( 1973) extended the 
workb v lZeFdovich & PikeFne r ( 1969) and lPenston & Browi] 
(Il970h to numerically compute isobaric flows i n thre e di- 
mensional spherical g eometry. Cowie & McKee il977h and 
iMcKee & Cowid ( 1197 7. hereafter MC77) considered a model 
of spherical clouds surrounded by the hot ionized medium 
and derived an analytic formula of the evaporation rate 
of clouds. MC77 also considered cold clouds surrounded 
by the WNM, which can be dir ectly compared with the 
work by 'Graham & Langei* ('1973V In our previous paper, 
Nagashima. Koyama & Inutsuka (2005), we considered the 
growth of spherical clouds and showed a new way to ob- 
taining approximate analytic formula of the evaporation and 
condensation rate. However, because we adopted a sim- 
ple polynomial form of the cooling function proposed by 
Elphick, Regev & Shaviv ( 1992) for simplicity, our under- 
standings were only qualitative. In this Letter, we shall show 
the evaporation rate for a realistic cooling function by using 
a full numerical simulation and compare it with the approxi- 
mate analytic formula and with the result of MC77. We would 
like to stress that our analysis can be applied to any spherical 
cold clouds surrounded by warm gas under pressure equilib- 
rium irrespective of their origin. Thus the results shown here 
are quite general. 

Our aim is thus to estimate the evaporation timescale of 
cold H I clouds as a function of the cloud size and the pressure 
of the ambient WNM. This paper is outlined as follows. In §2 
we compute the evaporation rate of cold clouds by using a full 
numerical simulation and compare it with MC77's. In §3 we 
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show the evaporation timescale. In §4 we discuss the depen- 
dence of the timescale on the size and pressure. We derive 
a new description of the evaporation rate as an improvement 
upon MC77. §5 we provide conclusions. 

2. EVAPORATION RATES 

Below we assume that clouds are spherically symmetric for 
simplicity, so that we compute only radial structure of clouds. 
Thus, the basic fluid equations are written as 
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r'^ dr ^ ' dr 
and the equation of state, p ~ pksT I p, where k{T) = 2.5 x 
erg ^ crn ^ K ^ s ^ is the conductivity for neutral 
gas ('Parker" 195 3"), ks the Boltzmann constant, p the mean 
molecular mass so that p = pn, and C is the heat-loss function 
defined as p£ = n^A — nV, where A and F are the cooling 
and heating rates, respe ctively. W e adopt a simple analytic 
fitting function given bv lKovama & Inutsuka (2002), 
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and r = 2 X 10 erg s ^, which are taking into account 
many processes, e.g., photoelectric heating from dust grains 
and PAHs, heating by cosmic rays and X- rays, and atomic line 
cooling from hydrogen Lya and C II (Ko vama & Inutsuka 
I2OOO). Throughout this paper, we ignore the effect of mag- 
netic fields. This might alter the results because the con- 
ductive heat flux is restricted along magnetic field lines, and 
therefore it should be taken into account in future analysis. 

Before performing full numerical simulations, we show 
the results for a quasi-steady state (QSS), under which we 
transform the time derivative to the spatial one, d/dt 
—Rcd/dr, where Rc is the cloud size and Rc is the velocity 
of the cloud size changing. Then the above equations become 
a set of ordinary differential equations. The details of the 
metho d to compute the structure will be foun d in a separate 
paper (INagashima. Kovama & Inutsukal2006l) . In Fig.l, sev- 
eral snapshots of evaporation rates, M(r) = Anr'^ p(r)v{r), 
are shown by thin dashed lines for the case of the saturation 
pressure, p = psat — 2823/0^. In contrast to the case of plane- 
parallel geometry, the cloud evaporates by the heat conduction 
and the size decreases as time passes. As we shall show in fj^ 
the motion of the front shown in this figure is purely driven by 
the curvature effect (Nagashima. Kovama & Inutsuka 2005). 

Next, we show the results obtained by full-numerical sim- 
ulations. _rhe_code_^s based on the second-order Godunov 
method (Ivan Leedll977h in the Lagrange coordinate. We im- 
pose a constant pressure at the outer boundary, that is, the 
WNM pressure is given as the boundary condition. The ini- 
tial structure of the cloud is computed by the QSS model with 
the size of about 0.1 pc. Correctly giving the initial condition 
is very important for this computation. If the deviation of the 
initial condition from the correct structure is large, spherical 
sound waves emerge and we could not extract the informa- 
tion on the evaporation. Several snapshots at which the size 
of the cloud becomes the same as that shown for the QSS are 
indicated by the thin solid lines in Fig.l. It is evident that the 
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Fig. 1 . — Evaporation rate, M for the case of the saturation pressure. Thin 
solid and dashed lines represent snapshots of M = 47rr^p(r)ti(r) as a func- 
tion of the distance from the cloud center for the full numerical simulation 
and the QSS calculation, respectively. Thick solid curve indicates M as a 
function of the cloud size. The size of clouds are defined in the text. Thick 
dashed curve indicates M given by MC77. 



results given by the full simulation agree well with those given 
by assuming the QSS probably because the Mach number of 
the flow is much below unity. Therefore we can consider that 
the QSS is a good assumption for the description of the evap- 
oration of clouds. 

The thick solid curve in Fig. 1 denotes the evaporation rate 
as a function of the cloud size, M{Rc), where the edge of 
the cloud is defined as a radius at which the density becomes 
("■CNM + "-wnm)/2. We have confirmed that the rate is pro- 
portional to the size, M{Rc) cx Rc- In fact, this is the same 
scaling against the size as that given by MC77. Their evapo- 
ration rate is 



Ml 



MC77 — 



lQ'KpK{Tf)Rc 



= 1.3 X IO'^'T]'^ ( — I gs 



(5) 

where T/ is the temperature of the WNM, Tj ~ 6.4 x 10^ K 
for the adopted heat-loss function. This rate is plotted by the 
thick dashed curve. When we write our evaporation rate in a 
similar way. 



M = 3.1 X lO^^T^^^^ 



gs 



(6) 



therefore this is about a factor of four smaller than MC77's 
in spite of the same dependence on the size. The form of the 
evaporation rate shall be discussed in fj^ in a different way 
from MC77. Note that the above rate is valid only for the case 
of the saturation pressure. In the next section, we show the 
evaporation timescale for other cases. 

3. EVAPORATION TIMESCALE 

Hereafter we define the evaporation timescale as icvap = 
M/ M, where the cloud mass M is estimated by integrating 
the mass density within the cloud radius, Rc, therefore Mc — 
47rpcNM-Rc/3. Fig. 2 shows the evaporation timescale against 
the cloud size Rc- The solid straight line indicates iovap for 
the case of the saturation pressure, p = psat, and for the QSS. 
The crosses on the solid line denote the result from the full 
numerical simulation. Similar to Fig.l, both are in excellent 
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Fig. 2. — Evaporation timescales defined as Mc/M as a function of tlie ra- 
dius of clouds. Solid, daslied, and dot-dashed lines indicate the timescale 
for p = psat, 2000fcB, and 40()0A:b, respectively, under the assumption 
of the QSS. Crosses denote simulation results. Note that Rc > 0.1 for 
Pwnm/^b = 4000 corresponds to the condensation timescale. 
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Fig. 3. — Critical radius i?crit plotted by the solid line against pressure of 
the WNM. Jeans length Rj is also plotted by the dashed line. 



agreement. It is easily confirmed that ^cvap oc i?^ oc M^/'^ 
forp = Psat, reflecting the previous result M cx Rc and M cx 
R^. Here it should be worth noting that tiny H I clouds with 
the size of ^ 10^^ pc must disappear in ^ IMyr irrespective 
of the ambient pressure. 

The dashed and dot-dashed lines indicate the evaporation 
timescale for the cases of p/fcs = 2000 and 4000 K cm^'^, 
respectively. Crosses on the lines also denote the results from 
the full numerical simulations. Again we can see the good 
agreement. Note that clouds always evaporate for p < psat, 
but large clouds can grow for p > psat, which is shown by 
the thin dot-dashed line at i? > 0.1 pc. This means that there 
exists a critical radius for growth (condensation), i?ciit (see 
the next section). Hereafter we use icvap also as the growth 
timescale. For these pressures, the size dependence of <cvap is 



cvap 

partly different from that for the saturation pressure. While it 
is very similar to that for the saturation pressure, icvap oc R^ 
at Rc <C i?crit, it becomes tovap Rc at Rc ^ -Rciit- This 
is very different feature from the expectation from MC77. In 
the next section, we explain what the size dependence comes 
from. 



4. DISCUSSION 

Below we derive an approximate evaporation rate in th e 
same way as that of Naga shima. Kovama & Inutsuka (2005!). 
Firstly, we assume that the evolution is isobaric because the 
motion of fluids is much slower than the sound speed. Then 
the energy equation Q can be described as an evolution equa- 
tion of enthalpy. 
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Next, again because of the slow motion of fluids, we assume 
that the fluids are in a QSS, so that the time derivative d/dt 
can be replaced with —Rcd/dr, 
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P"— = -P-C + -^K— + -K— , 
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— I j-i dr 



(8) 



where u= v—Rc- The first and second terms in the right hand 
side emerge independent of geometry of fronts. Therefore, if 
it is reasonable to assume that the structure of fronts is almost 
independent of geometry, we can replace the two terms with 
the left hand side for the case of the plane-parallel geometry, 

J kB dT 1 kB , .dT 2 dT 

7 — P"-^ - T — P^iiP)^ + -f^^: (9) 

7— 1/i or 7— 1/i or r or 

where ui{p) is the fluid velocity with respect to the front rest 
frame, and depends only on pressure. Noting the fact that the 
derivative of temperature, dT /dr has a non-zero, finite value 
only at the front, we can omit the derivative and substitute Rc 
for r, 

7 - 1 ^ 2 HiR 



u{R) ~ ui(p) 



(10) 



7 ks Rc PR 

where the subscript R stands for values at the front, r = Rc- 
To obtain a formula for the size evolution, we need a fur- 
ther assumption that a convergence factor from u to Rc is 
independent of geometry, / = —Rc/u. This becomes ex- 
act for the case of the thin-front limit. For three-dimensional 
spherical geometry, from the mass-flux conservation, we ob- 
tain — 47rr^pcNM^c = 4:TtRIpiiur, where v — inside the 
cloud. Then /sphere = (-RcPfl)/(''^PCNM). For the plane- 
parallel geometry, — pcnmc — Prur, where c is the speed of 
the front with respect to the cloud (CNM) rest frame, which 
corresponds to Rc- Then /p.p. — prJpcnm- If the front 
is sufficiently thin, and therefore the density p is pcNM even 
just inside the front, we obtain /sphere Pr/pcnm- This is 
equal to /p.p.. In reality, however, the front has a non-zero 
thickness. So the use of the same / is an approximation. 
Finally, we thus obtain the size evolution. 



Rr 



c{p) - / 



7^1^ 2 kr 
7 kB Rc PR ' 



(11) 



where we explicitly write the argument of c, p, to stress that 
it depends only on pressure. Forp > psat, c(p) > 0, and vice 
versa. This equation tells us the existence of a critical radius 
at which Rc — Q when c(p) > 0, 



Rn 



f 



j-l p 2 kr 



(12) 



7 kB c{p) PR ■ 
The evaporation rate is obtained from the above equation. 



M 



-inRlpcNuRc = Mp + Mc, 



(13) 
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where Mp and Afc are the "pressure" and "curvature" terms, 
respectively, and 



-47ri?^pcNMc(p) oc Rl, 



7 - 1 ^ 2 KR 

7 ks Rc PBa 



oc Rr 



(14) 
(15) 



From the different dependences on Rc, we can see that the 
pressure term dominates for large clouds, Rc ^ Rait, and 
that the curvature term does for small clouds, Rc <C i?crit- 
Using / ~ Pr/pcnm and 7 = 5/3, we obtain 



Mr. 



IGirpKuRc 



(16) 



5kB 

Thus the difference from MC77's result for the curvature term 



IS 



Mr 



MC77 




(17) 



because of k cx VT. Thus the difference between our result 
and MC77's can be partly explained by the above in addition 
to the thin-front approximation to derive the curvature term, 
Mr- MCll took into account only the curvature term. 

The size dependence of the evaporation timescale is thus 
dependent on the size. For small clouds, Rc ^ i?ciit, the 
timescale is icvap ~ M/ Mr oc R^. For large clouds, Rc ^ 
i?crit, it is tovap ~ M/Mp OC Rc- Thus the size dependence 
of the evaporation timescale shown in Figgis well explained 
by the above approximate evaporation rate we derived. 

A further important point is the relationship between the 
critical radius, i?crit, and the Jeans length for gravitational in- 
stability, i?j = csyT^jGp, where cs is the sound velocity. 
Fig|3l shows i?crit and i?j as a function of the ambient pres- 
sure. We find that there can exist clouds growing only by 
condensation for the case of p > psat- The timescale until 
the cloud size exceeds i?j is, however, quite long, as shown 
in Fig 121 

5. CONCLUSION 

We have investigated the evaporation rate and timescale of 
tiny H I clouds by using numerical simulations. We con- 
firmed that the results are almost recovered by assuming the 
evolution under the quasi-steady state because the fluid ve- 
locity is much slower than the sound velocity. We have found 



that clouds with a size of about 0.01 pc evaporate in approx- 
imately IMyr almost independent of pressure of the ambient 
WNM. This suggests that there might be some mechanisms 
to continuously form tiny H I clouds, or, that the ambient 
pressure around the clouds is much higher than the standard 
ISM pressure, if their existence is ubiquitous. The evapora- 
tion timescale of clouds larger than 0. 1 pc, however, depends 
strongly on the pressure. Clouds larger than a critical radius 
can even grow if the pressure is higher than the saturation 
pressure. 

In order to physically understand the simulation results, we 
derived an analytic formula for evaporation by assuming the 
isobaric evolution and a structure of the interface independent 
of geometry. The obtained evaporation rate has two separate 
terms: "pressure" and "curvature" terms. The former is inde- 
pendent of geometry, that is, it emerges even in analysis for 
plane-parallel geometry, and it becomes zero for the satura- 
tion pressure. The latter is proportional to the size of clouds, 
that is, the curvature. We have found that the obtained evap- 
oration rate is a natural extension of an evaporation rate ob- 
tained by MCll, which corresponds to our "curvature" term, 
but MC77's rate is a little higher than ours. For lower pressure 
than the saturation pressure, the signs of the two term are the 
same. Therefore clouds always evaporate. For higher pres- 
sure, on the other hand, the signs are different. Therefore a 
critical size exists and the size of clouds determines whether 
they grow or evaporate. We have confirmed that the simu- 
lation results show the same dependence on the size as the 
analytic formula. 

We would like to stress that our analysis can be adapted ir- 
respective of the origin of the tiny H I clouds. In order to ob- 
tain information on the formation process, statistical proper- 
ties such as mass function of clouds should be required. This 
will be done in future. 
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